Magnetization process in a chiral p-wave superconductor with multi-domains 
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A simulation study for the magnetization process is performed for the multi-domain state in a 
chiral p-wave superconductor, using the time-dependent Ginzburg-Landau theory. The external field 
penetrates inside as core-less vortices through the domain wall, forming the vortex sheet structure. 
We find that, with increasing magnetic fields, the domain walls move so that the unstable domains 
shrink to vanish. Therefore, the single domain structure is realized at higher fields. 
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A quasi-two dimensional superconductor Sr2Ru04 is 
an unconventional superconductor, and the pairing sym- 
metry is suggested to be a chiral p-wave pairing with 
the basic form p± — p x ± ip y and in-plane equal-spin 
pairing*^ The spin-triplet pairing is supported by the 
fact that the Knight shift does not change below the 
superconducting transition temperature T c . 3 Since the 
spontaneous moment appears below T c in the observation 
of the muon spin relaxation (/xSR), time-reversal symme- 
try of the pairing function is broken^ The internal field 
distribution of the vortex state, which is observed by a 
small angle neutron scattering (SANS), is consistent to 
the chiral p-wave pairing stated While the pairing func- 
tion of Sr2Ru04 may have an additional horizontal node, 
such as (p x ± ip y ) cosp z pi it is intrinsic in our study that 
the pairing function has a factor of the chiral component 
Px ± ipy 

In the chiral p-wave superconductor, a p+ state and 
a p_ state are degenerate in free energy. Therefore, the 
multi-domain structure may realize, i.e., some regions in 
a sample are p+ domains and others are p- domains. Be- 
tween the p+ and p- domains, domain walls appear as 
topological defects, which are not easily destroyed^ In 
the vortex state when a magnetic field H is applied to 
z direction, the degeneracy of the p+ and p_ states is 
removed. When H [[ z and charge e > (or equivalently 
when H || — £ and e < 0), the free energy of the p- state 
is lower than that of the p+ stated because the vortex 
structures are different between the p + and p_ domains. 
In the p- (p+) state, the opposite p + (p_) component 
induced around the vortex core has large (small) ampli- 
tudeAifliii Since the p- state is stable in the vortex state, 
it is interesting to see how the multi-domain structure at 
a zero field changes to the single p_ domain, when a mag- 
netic field is applied. We note that the chirality depen- 
dence is defined relative to the magnetic field direction, 
the free energy of the p + state becomes smaller when 
the magnetic field is applied to the reverse direction —z. 
Study of vortices trapped at the domain wall is also im- 
portant since they are considered to have a strange struc- 
ture called a "vortex sheet" ji£ where half flux-quantum 
vortices are aligned along the domain wall&iiiiii^i& 

The purpose of this paper is to study the magneti- 
zation process in order to understand how the multi- 
domain structure changes to the single p_ domain by 



applying magnetic fields in the chiral p-wave supercon- 
ductor. We also investigate the roles and properties of 
the vortex sheet structure at the domain wall appearing 
in the magnetization process. The static properties of the 
vortex or domain wall structure were studied also by the 
two-component Ginzburg-Landau (GL) theoryi^&iSiiLiS 
For the study of the magnetization process, we use the 
time-dependent Ginzburg-Landau (TDGL) theory as a 
phenomenological approach in our qualitative study, with 
expectation that vortices move so as to approach the 
free energy minimum state. The magnetization process 
of a conventional superconductor was studied by Kato 
et alM^ The present authors studied the static and 
dynamical properties of the vortex sheet structure in 
a two-component superconductor based on the TDGL 
theoryiiii4 In this paper, we use the GL free energy for 
the chiral p-wave superconductor, and study the magne- 
tization process. 

To obtain the two-component GL equation in the chiral 
p-wave superconductor, the pair potential is decomposed 
as A(r, p) = ?7 + (r)0 + (p) + r/_(r)^_(p) with the order 
parameter r]±(r), where r is the center-of-mass coordi- 
nate of the Cooper pair. The pairing functions <fr±(p), 
depending on the relative momentum p of the pair, are 
given by the chiral p-wave type such as p x ± ip y . The GL 
free energy density is written as 
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+2^+1 V| 2 + d + V * + 2 V 2 _ ) + V * + (q 2 x + q 2 y ) V+ 

+c 3 (vUlv- + v-q 2 -v+) (i) 

in the dimensionless form^ where q± = q x ± iq y , q = 
(Ti/i)V — 27r(2e//ic)A with vector potential A. hc/2\e\ = 
(f>o is a flux-quantum. The coefficients are related to the 
pairing function and the Fermi surface structure as 
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where v± = (v x ± iv y )/2 with a Fermi velocity (v x ,v y ), 
and (• ■ ■} indicates the average on p along the Fermi sur- 
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face. C\ and C3 are anisotropy parameters which are 
finite when the pairing functions or the Fermi surface 
have fourfold symmetric structure. For the isotropic case, 
C\ = C3 = due to the vanishing Fermi surface average 
in Eq. (J5J. As the detailed forms of <j>± and the Fermi 
surface structure have not been established yet, we treat 
the coefficients in Eq. @ as arbitrary parameters. In 
this study we set C\ = C3 = for simplicity to exclude 
the additional anisotropy effect. We show results calcu- 
lated for C 2 = 0.3 and T = 0.5T C . 

In our simulations, we use the TDGL equation coupled 
with Maxwell equationSfi2i21i 
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B = V x A. 
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The supercurrent j s = (j s , x ,j s , y ) oc (df /dA x ,df /dA y ) 
is given by j s . x = Re[r)+{q x r] + ) + r]*_(q x r)^) + 

C2{vUv-v-)+v-{q+v+)}+C3{v*+{q+v-)+v-{q-v+)}}, 
kv = Re [vl(q v v+) + v-(qyv-) - i^iv+ii-v-) - 

V-(<1+ 7 1+)} + iCsiv+il+V-) - V-(<l-V+)}}- The length, 
field, and time are, respectively, scaled by the coherence 
length £ , H c2 .o = 4>q/2itQ, and t = 47r£ 2 K 2 er/c 2 with 
the normal state conductivity cr^SS However, we here 
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scale rj± by 770 instead of t]q(T) = 770(1 — T/T c ) 
a uniform solution of 774- when 77 _ = and T = 0. The 
calculations are performed in a two-dimensional rectan- 
gular area with a size 200£o x 100£o- Outside the open 
boundary, we set 77+ = r/_ = and B(r) — H with an 
applied field H. We set the GL parameter k = 2.7. 

Our calculation of the magnetization process for the 
multi-domain state is shown in Fig. ^1 At a zero field, 
we prepare the state where the right-hand side half re- 
gion is a p + state (|?7+| ~ 1, 77- =0), and the left-hand 
side half region is a p_ state (|r?-| ~ 1, 77+ = 0). A 
straight domain wall appears at the center between the 
p + and p_ domains. The relative phase of 77+ and ?7_ 
is 7r, which minimizes the free energy in our case. We 
increase H gradually with a slow rate SH/St = 5 x 10~ 6 . 
The left panels in Fig. ^ show the color-density plot of 
|?7+(r)| and |r7_(r)|. A green (red) region indicates the 
p + (p-) domain. The center panels show ^(r)) and 
1 77^, (r) I , when we define the p x component 77^ and the p y 
component 77^ as A(r,p) = rj x (r)(f> x (p) + ri y (r)<f> y (p) with 
4>x = (4>++<t>-)/2 ~p x and^y = ((/j + -</3_)/2i ~ p y . Yel- 
low region indicates that 1 77^ (r) | ~ |%(r)| ~ 1 in the p + 
or p_ domains. The right panels represent the internal 
field distribution B(r), which can be observed directly. 

At low fields in the Meissner state [Fig. Ufa)], magnetic 
fields penetrate inside through the domain wall, forming 
the vortex sheet structure. A straight domain wall at 
H = begins to meander by the penetration of vortices 
at finite fields. To see the vortex sheet structure, we show 
the spatial structure in Fig. |3 magnifying the enclosed 

area in Fig. Qfa). In the amplitude 1 77 (i*) | and |?7+( r )l m 

Figs. |U^a) and[2{b), we do not see the singularity of the 



vortex center. It is because the singularity points with 
the phase winding 2-7T of vortices are at the opposite re- 
gion across the domain wall, where the amplitude of the 
order parameter is well suppressed. That is, the winding 
center of 77+ (?7— ) is located in the p- (p+) domain, as 
shown by a solid circle in Fig. [2Kb) (Fig. Ufa)). There- 
fore, vortices at the domain wall are core-less vortices. 

When we consider 77^ and 77^ instead of 77+ and 77- , as 
shown in the center panels in Fig. ^ the core-less vortices 
of the vortex sheet are seen as different structures. At 
a zero field, the domain wall is presented as a green line 
in the color density plot, since rj y ~ 1 and 77^ ~ due 
to the sign change of 77^ at the domain wall. When the 
magnetic field penetrates, as shown in the center panel 
of Fig. ^a) or in Figs. \%Lc) and^d), vortices of r/ y (red 
circle) enter along the domain wall from the boundary. 
These r\ y vortices are located slightly at the p+ domain 
side in the domain wall region. The green line between 
the 77y vortices changes to the vortex of 77^ (green circle) , 
when an inter- vortices distance of the 77^ vortices becomes 
short with increasing H . That is, the order parameters 77^ 
and T] y have vortex cores with a winding 2tt at different 
positions. This vortex sheet structure with the 77 x and 
r] y vortices alternately aligning along the domain wall is 
the same vortex sheet structure reported in our previous 
workii^ii^ The B(r) distribution has a ridge along the 
domain wall, and has a peak at the r/ y vortices, as shown 
in the right panel of Fig. ^a) or in Fig. Hfe). 

With increasing H , first, vortices penetrate inside only 
in the p_ domain [Fig. IHb)]. Later at higher field, vor- 
tices penetrate also in the p+ domain [Fig. [He)] . These 
indicate that the lower critical field H c \ in the p- domain 
is lower than that in the p + domain. This corresponds to 
the fact that the upper critical field H C 2 in the p- state is 
higher than that of the p+ stated The amplitude of the 
induced p+ component around vortices in the p- state is 
larger than that of the induced p_ component in the p+ 
state. 9 Since H c \ is related to the creation energy of the 
vortex core, therefore, the p- state with smaller creation 
energy of the vortex core has smaller H c i compared to the 
p+ stated In Figs. |Jb) and^c), we see vortices at the 
boundary region. This is because the penetrating mag- 
netic field decreases towards inside of the superconductor 
in the length scale of the penetration depth. Near H c i 
when H increases, vortices first appear at the boundary 
region where locally B(r) > H c \. 

With increasing the number of vortices along the do- 
main wall, the domain wall line moves meanderingly. 
With further increasing H, the domain wall moves so 
that the area of the p + domain shrinks as is seen in Fig. 
Did), where the vortex sheet structure still appears along 
the meandering domain wall. Finally at enough high 
fields, the p+ domain vanishes, and the single p- domain 
state is realized. 

Around the vortex cores in the p+ domain and the 
p_ domain, the opposite chiral component is induced as 
discussed in previous works AiSiii We also observe some 
double-winding An vortex near the boundary regions in 
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FIG. 1: (Color) Magnetization process for the multi-domain state in a 200£o x 100£o area with the open boundary. We start 
from the zero-field state where the left (right) hand side half region is a p+ (p~) domain. The applied field H at the boundary 
is gradually increased, (a) H/H C 2.o = 0.08, (b) 0.11, (c) 0.13 and (d) 0.32. The left panels show the color density plots of 
|»7+(r)| and |jj_(r)|. The center panels are for |r/a;(r)| and |%(r)|. The right panels present the internal field distribution B(r). 



the p- domain. They are seen as large green circles in 
the red region in the color density plots of \r] + \ and |t/_| 
in Figs. db)-^d), since the vortex core of the double- 
winding vortex in the ?y_ component is filled by the in- 
duced r] + component with a winding 0. The population of 
the double- winding vortex increases with increasing C2, 
which is the strength of the gradient coupling between 
the p+ and p_ components. In the B(r) distribution, 
the double-winding vortex has a local minimum of B(r) 
at the vortex center within the enhanced B(r) region 
around the vortex core, due to the induced 77+ compo- 
nent. Because this B(r) structure is not consistent to 
that observed by SANSf^ however, most of the vortices 
in Sr2Ru04 are not double-winding vortices. 

We also examine the single domain case of a p+ state 
at H = 0. Since the p^ state has lower free energy in 
the vortex state, the p+ domain has to change to the p_ 
domain by applying fields. With increasing H, vortices 
penetrate inside in the single p+ domain at H ~ 0.11. 
At a higher field H ~ 0.25, some small p_ domains are 
created at the boundary. With further increasing H, the 
area of the p- domain extends inside, as shown in Fig. 



13 Finally, the p+ domain shrinks to vanish also in this 
case, realizing the single p_ domain state. 

Lastly, we compare the magnetization curves in Fig. 0] 
for three cases; (i) the multi-domain case when the p+ 
domain and the p_ domain coexist at low fields as shown 
in Fig. 2] (ii) the case of the single p+ domain at low 
fields as shown in Fig. |3 and (iii) the case of the single 
p_ domain. In the single domain case, B = in the 
Meissner state at H < H c i, and B appears in the mixed 
state at H > H c i due to the penetration of vortices. We 
find that H c \(j>- state) < H c i(p + state) also in Fig. 0] 
In the multi-domain case (thin solid line in Fig. 0, B is 
small but finite even in the Meissner state, because the 
magnetic field penetrates along the domain wall. In this 
case, the slope of B curve changes both at H c i(p + state) 
and at H c i(p_ state). 

In summary, we performed the simulation of the mag- 
netization process in a chiral p-wave superconductor, us- 
ing the TDGL theory with two components of the p+ 
and p- states. In the multi-domain case of the p+ and 
p_ domains, the magnetic field penetrates as core-less 
vortices along the domain wall, forming the vortex sheet 
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FIG. 2: Vortex sheet structure at H/H C 2,o = 0.08 in the 
enclosed area in Fig. 1(a). (a) |r?-(r)|, (b) |»7+(r)[, (c) |r?a;(r)|, 
(d) |%(r)|, and (e) B(r) are presented. Solid circles in (a) 
and (b) indicate the position of a phase winding 2n for r\- (r) 
and rj + (r), respectively. 




FIG. 3: (Color) Magnetization process in the single domain 
case of a p+ state at H — 0. Color density plots of |7j+(r)| and 
| J7_ (r) | are presented at (a) H/H C 2,o = 0.31, and (b) 0.35. 



external fields, the domain wall meanderingly moves, so 
that the area of the p+ domain shrinks. Then, the un- 
stable p + domain vanishes at a high field, and the single 
domain of the stable p- state is realized. Even in the 
case of the single p+ domain at a zero field, the p+ do- 
main changes to the single p_ domain at a high field as 
in a similar way, after small p- domains are created at 
the boundary. 

Recently, anomalous internal field distribution due to 
the vortex coalescence was observed i2I These may be re- 
lated to the intrinsic character of the chiral p-wave su- 
perconductor, such as a domain structure. 
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FIG. 4: Magnetization curve at low fields. Averaged flux 
density B is plotted as a function of an external field H. We 
normalize H and B by H C 2,o- When the spatial average B 
is calculated from B(r), we exclude the boundary region of 
the width 10£o not to consider the surface effect at the sample 
boundary. A thin solid line is for the multi-domain case shown 
in Fig. 1. A thick solid (dashed) line is for the single domain 
case when we start from a single p+ (p~) domain at H — 0. 
The dash-dotted line shows B = H. 



structure, even in the Meissner state. With increasing 
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